It looks like the ECS distribution is not normally distributed, even though $B$ is.
👉 How does $\overline{\text{ECS}(B)}$ compare to $\text{ECS}(\overline{B})$? What is the probability that $\text{ECS}(B)$ lies above $\text{ECS}(\overline{B})$?
(MIT students only)
Please see the link for hw 9 transcript document on Canvas. We want each of you to correct about 500 lines, but don’t spend more than 20 minutes on it. See the the beginning of the document for more instructions. :point_right: Please mention the name of the video(s) and the line ranges you edited:
In the lecture notebook we introduced a mutable struct EBM (energy balance model), which contains:
the parameters of our climate simulation (C, a, A, B, CO2_PI, α, S, see details below)
a function CO2, which maps a time t to the concentrations at that year. For example, we use the function t -> 280 to simulate a model with concentrations fixed at 280 ppm.
EBM also contains the simulation results, in two arrays:
T is the array of tempartures (°C, Float64).
t is the array of timestamps (years, Float64), of the same size as T.
Hint
@bind log_CO2 Slider(❓)CO2 = 10^log_CO2We talked about two emissions scenarios: RCP2.6 (strong mitigation - controlled CO2 concentrations) and RCP8.5 (no mitigation - high CO2 concentrations). These are given by the following functions:
18.S191, fall 2020
👉 Find the lowest CO₂ concentration necessary to melt the Snowball, programatically.
homework 9, version 0
We talked about a second theory – a large increase in CO₂ (by volcanoes) could have caused a strong enough greenhouse effect to melt the Snowball. If we imagine that the CO₂ then decreased (e.g. by getting sequestered by the now liquid ocean), we might be able to explain how we transitioned from a hostile Snowball Earth to today's habitable "Waterball" Earth.
In this exercise, you will estimate how much CO₂ would be needed to melt the Snowball and visualize a possible trajectory for Earth's climate over the past 700 million years by making an interactive bifurcation diagram.
In the lecture notebook (video above), we had a bifurcation diagram of $S$ (solar insolation) vs $T$ (temperature). We increased $S$, watched our point move right in the diagram until we found the tipping point. This time we will do the same, but we vary the CO₂ concentration, and keep $S$ fixed at its default (present day) value.
Hint
The function findfirst might be helpful.
Submission by: Jazzy Doe (jazz@mit.edu)
Abstraction, lines 1-219; Array Basics, lines 1-137; Course Intro, lines 1-144 (for example)
Answer:
Hello world!
Hint
Use a condition on the albedo or temperature to check whether the Snowball has melted.
If you like, make the visualization more informative! Like in the lecture notebook, you could add a trail behind the black dot, or you could plot the stable and unstable branches. It's up to you!
Reveal answer:
Again, look inside simulated_model and notice that T and t have accumulated the simulation results.
In this simulation, we used T0 = 14 and CO2 = t -> 280, which is why T is constant during our simulation. These parameters are the default, pre-industrial values, and our model is based on this equilibrium.
👉 Run a simulation with policy scenario RCP8.5, and plot the computed temperature graph. What is the global temperature at 2100?
Additional parameters can be set using keyword arguments. For example:
Model.EBM(14, 1850, 1, t -> 280.0; B=-2.0)Creates the same model as before, but with B = -2.0.
👉 Create a slider for CO2 between CO2min and CO2max. Just like the horizontal axis of our plot, we want the slider to be logarithmic.
You can set up an instance of EBM like so:
👉 In what year are we expected to have doubled the CO₂ concentration, under policy scenario RCP8.5?
In lecture 21 (see below), we discovered that increases in the brightness of the Sun are not sufficient to explain how Snowball Earth eventually melted.
Have look inside this object. We see that T and t are initialized to a 1-element array.
Let's run our model:
We are interested in how the uncertainty in our input $B$ (the climate feedback paramter) propagates through our model to determine the uncertainty in our output $T(t)$, for a given emissions scenario. The goal of this exercise is to answer the following by using Monte Carlo Simulation for uncertainty propagation:
mimetext/htmlrootassigneelast_run_timestampA <persist_js_state·has_pluto_hook_features§cell_id$06c5139e-252d-11eb-2645-8b324b24c405depends_on_disabled_cells§runtime bpublished_object_keysdepends_on_skipped_cells§errored$b9f882d8-266b-11eb-2998-75d6539088c7queued¤logsrunning¦outputbodymimetext/plainrootassigneelast_run_timestampA <$ypersist_js_state·has_pluto_hook_features§cell_id$b9f882d8-266b-11eb-2998-75d6539088c7depends_on_disabled_cells§runtime&published_object_keysdepends_on_skipped_cells§errored$1843dae6-2689-11eb-2aaf-036eb2f9341dqueued¤logsrunning¦outputbodyprefixFloat64elements 2.08467text/plain5.0928text/plain2.38992text/plain4.42845text/plain4.36886text/plain2.22342text/plain1.53348text/plain3.221text/plain 2.86569text/plain 2.37656text/plain 1.6003text/plain 2.44421text/plain 1.94336text/plain2.47175text/plain2.75275text/plain4.62015text/plain2.13501text/plain1.56871text/plain2.62535text/plain6.55421text/plainmore~2.73358text/plain2.16text/plain2.92945text/plain3.338text/plain2.98041text/plain1.88202text/plain2.5971text/plain2.34947text/plain2.13581text/plain2.82829text/plaintypeArrayprefix_shortobjectid805657125eb905a0mime!application/vnd.pluto.tree+objectrootassigneeECS_sampleslast_run_timestampA E]persist_js_state·has_pluto_hook_features§cell_id$1843dae6-2689-11eb-2aaf-036eb2f9341ddepends_on_disabled_cells§runtime|ppublished_object_keysdepends_on_skipped_cells§errored$35f87c2e-253d-11eb-0d79-61d89c1d9b5equeued¤logsrunning¦outputbody281.8382931264455mimetext/plainrootassigneeCO2last_run_timestampA Fspersist_js_state·has_pluto_hook_features§cell_id$35f87c2e-253d-11eb-0d79-61d89c1d9b5edepends_on_disabled_cells§runtime;published_object_keysdepends_on_skipped_cells§errored$e296c6e8-259c-11eb-1385-53f757f4d585queued¤logsrunning¦outputbody1👉 What is the probability that we see more than 2°C of warming by 2100 under the low-emissions scenario RCP2.6? What about under the high-emissions scenario RCP8.5?
👉 Change the value of $B$ using the slider above. What does it mean for a climate system to have a more negative value of $B$? Explain why we call $B$ the climate feedback parameter.
Before you submit
Remember to fill in your name and Kerberos ID at the top of this notebook.
A recent ground-breaking review paper produced the most comprehensive and up-to-date estimate of the climate feedback parameter, which they find to be
$$B \approx \mathcal{N}(-1.3, 0.4),$$
i.e. our knowledge of the real value is normally distributed with a mean value $\overline{B} = -1.3$ W/m²/K and a standard deviation $\sigma = 0.4$ W/m²/K. These values are not very intuitive, so let us convert them into more policy-relevant numbers.
Definition: Equilibrium climate sensitivity (ECS) is defined as the amount of warming $\Delta T$ caused by a doubling of CO₂ (e.g. from the pre-industrial value 280 ppm to 560 ppm), at equilibrium.
At equilibrium, the energy balance model equation is:
$$0 = \frac{S(1 - α)}{4} - (A - BT_{eq}) + a \ln\left( \frac{2\;\text{CO}₂_{\text{PI}}}{\text{CO}₂_{\text{PI}}} \right)$$
From this, we subtract the preindustrial energy balance, which is given by:
$$0 = \frac{S(1-α)}{4} - (A - BT_{0}),$$
The result of this subtraction, after rearranging, is our definition of $\text{ECS}$:
$$\text{ECS} \equiv T_{eq} - T_{0} = -\frac{a\ln(2)}{B}$$
Properties of an EBM obect:
| Name | Description |
|---|---|
A | Linearized outgoing thermal radiation: offset |
B | Linearized outgoing thermal radiation: slope. or: climate feedback parameter |
α | Planet albedo, 0.0-1.0 |
S | Solar insulation |
C | Atmosphere and upper-ocean heat capacity |
a | CO₂ forcing effect |
CO2_PI | Pre-industrial CO₂ concentration |
Before working on the homework, make sure that you have watched the first lecture on climate modeling 👆. We have included the important functions from this lecture notebook in the next cell. Feel free to have a look!
👉 Write a function step_model! that takes an existing ebm and new_CO2, which performs a step of our interactive process:
Reset the model by setting the ebm.t and ebm.T arrays to a single element. Which value?
Assign a new function to ebm.CO2. What function?
Run the model.
👉 Generate a probability distribution for the ECS based on the probability distribution function for $B$ above. Plot a histogram.
The albedo feedback is implemented by the methods below:
mimetext/htmlrootassigneelast_run_timestampA👉 Write a function temperature_response that takes a function CO2 and an optional value B as parameters, and returns the temperature at 2100 according to our model.
We used two helper functions:
👉 Create a graph to visualize ECS as a function of B.
mimetext/htmlrootassigneelast_run_timestampA <persist_js_state·has_pluto_hook_features§cell_id$aed8f00e-266b-11eb-156d-8bb09de0dc2bdepends_on_disabled_cells§runtime published_object_keysdepends_on_skipped_cells§errored$269200ec-259f-11eb-353b-0b73523ef71aqueued¤logsrunning¦outputbody!To compute ECS, we doubled the CO₂ in our atmosphere. This factor 2 is not entirely arbitrary: without substantial effort to reduce CO₂ emissions, we are expected to at least double the CO₂ in our atmosphere by 2100.
Right now, our CO₂ concentration is 415 ppm – 1.482 times the pre-industrial value of 280 ppm from 1850.
The CO₂ concentrations in the future depend on human action. There are several models for future concentrations, which are formed by assuming different policy scenarios. A baseline model is RCP8.5 - a "worst-case" high-emissions scenario. In our notebook, this model is given as a function of $t$.
👉 Inside the plot cell, call the function step_model!.
Hello world!
Fantastic!
Splendid!
Great!
Yay ❤
Great! 🎉
Well done!
Keep it up!
Good job!
Awesome!
You got the right answer!
Let's move on to the next section.
$B =$
Hello world!
Below we have an empty diagram, which is already set up with a CO₂ vs $T$ diagram, with a logirthmic horizontal axis. Now it's your turn! We have written some pointers below to help you, but feel free to do it your own way.
Let's create a package environment:
Hint
Start by writing a function equilibrium_temperature(CO2) which creates a new EBM at the Snowball Earth temperature T = -48 and returns the final temperature for a given CO2 level.
👉 Does accounting for uncertainty in feedbacks make our expectation of global warming better (less implied warming) or worse (more implied warming)?
Just some helper functions used in the notebook.
👉 What happens when $B$ is greater than or equal to zero?
The climate feedback parameter $B$ is not something that we can control– it is an emergent property of the global climate system. Unfortunately, $B$ is also difficult to quantify empirically (the relevant processes are difficult or impossible to observe directly), so there remains uncertainty as to its exact value.
A value of $B$ close to zero means that an increase in CO₂ concentrations will have a larger impact on global warming, and that more action is needed to stay below a maximum temperature. In answering such policy-related question, we need to take the uncertainty in $B$ into account. In this exercise, we will do so using a Monte Carlo simulation: we generate a sample of values for $B$, and use these values in our analysis.
mimetext/htmlrootassigneelast_run_timestampAThe plot below provides an example of an "abrupt 2xCO₂" experiment, a classic experimental treatment method in climate modelling which is used in practice to estimate ECS for a particular model. (Note: in complicated climate models the values of the parameters $a$ and $B$ are not specified a priori, but emerge as outputs of the simulation.)
The simulation begins at the preindustrial equilibrium, i.e. a temperature $T_{0} = 14$°C is in balance with the pre-industrial CO₂ concentration of 280 ppm until CO₂ is abruptly doubled from 280 ppm to 560 ppm. The climate responds by warming rapidly, and after a few hundred years approaches the equilibrium climate sensitivity value, by definition.
Properties of an EBM obect:
| Name | Description |
|---|---|
A | Linearized outgoing thermal radiation: offset |
B | Linearized outgoing thermal radiation: slope. or: climate feedback parameter |
α | Planet albedo, 0.0-1.0 |
S | Solar insulation |
C | Atmosphere and upper-ocean heat capacity |
a | CO₂ forcing effect |
CO2_PI | Pre-industrial CO₂ concentration |